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Standard Monte Carlo simulations have been performed on improved lattices to determine the 
wave functions and the sizes of the scalar and tensor glueballs at four lattice spacings in the range 
a = 0.05 — 0.145 fm. Systematic errors introduced by the discretization and the finite volume are 
studied. Our results in the continuum limit show that the tensor glueball is approximately two 
times as large as the scalar glueball. 



PACS numbers: 



I. INTRODUCTION 



Lattice calculations have yielded very accurate spectro- 
scopic information regarding low-lying hadrons 0,12 and 
glueballs [E EL IE IE Ll HI ■ These calculations show that 
the lowest-lying scalar, tensor and axial vector glueballs 
lie in the mass region 1-2.5 GeV. While there is a long 
history of hadron and glueball mass calculations in lat- 
tice QCD, progress in determining their wave functions 
has not been so rapid. Although the spectrum of pure 
gauge QCD cannot be directly compared to experiment, 
it is important to pin it down because so little is known 
about glueball properties. The mixing of glueballs with 
qq mesons due to the presence of light dynamical quarks 
makes it difficult in distinguishing them from ordinary 
mesons, both in full lattice QCD and in experiments 
Enormous effort have been made in pursuit of glueballs 
and their properties ^3] • Accurate lattice calculations of 
their size, matrix elements and form factors would help 
considerably in the attempt to realize their experimental 
identification. 

Glueball wave functions and sizes have been studied 
previously [U 02 d , but much of the early works 
contain uncontrolled systematic errors, most notably 
from discretization effects. However, calculations using 
operator overlaps obtained from variational optimization 
for an improved lattice gauge action showed that the 
scalar and tensor glueballs are of typical hadronic dimen- 
sions. A straightforward procedure to determine glueball 
size is to compute measure the glueball wave function, 
much in the same way as in the meson and baryon wave 
functions were measured. 

Although Ref. produced interesting results, the 

approach used there suffers from a basic problem: The 
observables are calculated from lattice version of the 2- 
glue operator, which risks the mixture of glueball states 
with the flux state 1 . Also the results arc of limited inter- 



est, because of their manifest dependence on the gauge 
chosen. In this study we take a more direct approach to 
the problem: Instead of fixing a gauge or a path for the 
gluons, we measure the correlation functions from spa- 
tially connected Wilson loops which, being the expecta- 
tion values of closed-loop paths, are gauge invariant. The 
procedure has been used in the measurement of the wave 
functions of the scalar and tensor glueballs in the lattice 
Monte Carlo calculations in 3-dimensional U(l) theory 
[l5|. This approach has the advantage that the detan- 
gling of the glueball and torelon is usually taken care of 
automatically by the choice of Wilson or Polyakov loops. 
This means that glueball operators cannot create a sin- 
gle torelon state, but the creation of torelon-antitorelon 
pairs is possible. 



In this paper, we demonstrate the efficiency of our 
method in calculating SU(3) glueball wave functions us- 
ing the tadpole-improved Symanzik action 6]. Since we 
want to explore the nature of the wave function, we fo- 
cus on the sizes and masses of two of the lighter SU(3) 
glueball states, the scalar and the tensor. To avoid mix- 
ing effects, we adopt a quenched lattice QCD as a nec- 
essary first step before attempting to include the effect 
of dynamical quarks. It is worth mentioning here that 
quenched lattice QCD accurately predicts the impor- 
tant nonperturbative quantities as well as the masses of 
hadrons, mesons, and baryons. The remaining content 
of this paper are organized as follows. In Sec. [n]we give 
a prescription for calculating the glueball wave function 
and size based on the information from glueball correla- 
tors with the smearing method. We present and discuss 
our results in Sec. IIIII We close section with a compari- 
son of our results with related QCD results. Concluding 
remarks are given in Sec. IIVI 



1 The link-link operator used there sums up a large number of 



loops; some of these loops have a zero winding numbers and 
project onto glueballs others have non-zero winding numbers and 
project onto flux states, also called torelons. 
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II. WAVE FUNCTIONS OF GLUEBALLS 

The lattice observables are computed as follows. First, 
we calculate the lattice operator 

$,•(?,*) = 5^ + (1) 

X 

where 4> is the plaquette operator and $ is the two- 
plaquette or two-loop component of the glueball wave 
function. The summation over x is performed in order 
to obtain the zero-momentum projection. The r depen- 
dence is reflected in the length of the links required to 
close the loops. 

Given a pair of orbit and spin quantum numbers L 
and S for the desired glueball, one can construct oper- 
ators W(\ f \,t) with J = 0,1,2, P = ±, and C = ± 
from suitable linear combinations of the rotation, parity 
inversions and real and imaginary parts of the operators 
involved in $. In this study, we only consider S = L = 
(scalar (0++) and S = 2, L = (tensor 2++) glueball 
states. The two lattice observables measured are there- 
fore 

12 + $13 + $23] j (2) 

r 

W 2 (\f\,t) = ^ [2$ 12 - $ 23 - $13] . (3) 

r 

Using the glueball operators W, we consider the cor- 
relator 

C(r,t) = (Wi(r,t)W(r,0)), (4) 

where the vacuum expectation value must to be sub- 
tracted for the ++ glueball. The source can be held 
fixed, while the sink takes on the r dependence. This 
proves to be helpful in maintaining a good signal. In 
general, to measure the ground-state mass from the cor- 
relator, one seeks the region where the overlap with the 
lowest state is maximum and contributions from excited 
states are negligible. In principle, such a region always 
exists for large Euclidean time t in the zero-temperature 
case. However, in practice, it is difficult to use such large 
t in lattice QCD Monte Carlo calculations, since the cor- 
relator decreases exponentially with t and becomes so 
small for large t that it is comparable to its statistical 
errors. Hence, for the ground-state mass measurement, 
it is important that the ground-state overlap be suffi- 
ciently large. This leads to what amounts to a discrete 
search among wave functions. In the case of glueballs in 
quenched SU(3) lattice QCD, the overlap of the opera- 
tors given in Eqs. (2) and (3) is quite small as long as 
they are constructed from the simple plaquette operators. 
This small overlap originates from the fact that the "size" 
of the plaquette operator is smaller (specifically of 0{a)) 
than the physical size of the glueball ^(|. To improve 
the glueball operator, to enhance the ground-state con- 
tribution, we exploit the APE link smearing techniques 



|16lll7j to generate extended operators with sizes that are 
approximately the same as the physical size of the glue- 
ball. The smearing procedure is implemented through 
the iterative replacement of the original spatial link vari- 
able by a smeared link. For the smearing parameters, 
which play an important role in extracting the ground- 
state contribution, we used a smearing fraction of a = 0.5 
and a smearing number 7V smr = 15 in the present calcula- 
tion. The optimum smearing is determined by examining 
the ratio 

C(r,t + 1)/C(r, t) 

which should be close to unity for good ground sate dom- 
inance. In fact, in the range 4 < N sml < 11, the above 
ratio is almost unchanged. For our measurements, with 
a = 0.5 fixed, a typical value that proved to be sufficient 
for most cases was N smx = 4 (although the results ob- 
tained with this value are almost the same as those with 
iV sm r = 11). For fixed a, N smI plays the role of extend- 
ing the size of the smeared operator and hence can give 
the rough estimate of the physical size of the glueball. 

A second pass was made to compute the optimized 
correlation matrices 

C«(t) = (W(n,t)W( rj ,0)) - (W( ri ))(W( rj )). (5) 

Let ip^ be the radial wave function of the fc-th eigenstate 
of the transfer matrix. Then we have 

C y (t) =J2 a k<P (k) (n)i> {k) (rj)e- mkt - (6) 

k 

The glueball masses and the wave functions are extracted 
from the Monte Carlo average of Cy (t) by diagonalizing 
the correlation matrices C{t) for successive times t, 

C{t) = R(t)D(t)R(t), (7) 

where D is a diagonal matrix whose elements are eigen- 
values and R is a rotation matrix whose columns are the 
eigenvectors of C. Each eigenvector of C corresponds to 
an eigenstate (r) of the complete transfer matrix. 

In order to investigate the size of the glueball, the 
Bethe-Salpeter amplitude provides a convenient tool. In 
the nonrelativistic limit, it is expected to reduce to 
the glueball wave function in the first-quantized picture. 
This propert y c an be exploited to estimate the size of 
the glueball [Ij. Similar to the case for mesons, the 
wave function is expected to decrease exponentially with 
the separation r and is therefore fitted with the simple 
form 

V>(r) ee e - r/ro (8) 

to determine the effective radius r<j. 

The signal in the connected 2-point correlator C (t) for 
all observables does not last long to make to an expo- 
nential fit. For this reason, the data are presented in 
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terms of the effective mass read directly from the largest 
eigenvalue corresponding to the lowest energy, 



^ff 



Ao(r = 0,t=l) 



A (r = 0,t = 2) 



(9) 



However, to ensure the validity of our results, we com- 
pared them to those obtained using 



id 



off : 



Ao(t-l)-A (t) 
A («)-Ao(t + l) 



(10) 



III. 



SIMULATIONS RESULTS AND 
DISCUSSION 




An ensemble of gauge configurations were generated 
using a combination of the Cabibbo-Marinari (CM) algo- 
rithm and the overrelaxed method. Configurations were 
given a hot start and then 100 compound sweeps (we de- 
fine one compound sweep as one CM update followed by 
five over-relaxation sweeps) in order to equilibrate. Af- 
ter thermalization, configurations were stored every 50 
compound sweeps for 750 configurations. Measurements 
made on the stored configurations were binned into 10 
blocks with each block containing an average of 75 mea- 
surements. The mean and the final errors were obtained 
using the single-elimination jackknife method with each 
bin regarded as an independent data point. Four sets 
of measurements were taken to check the scaling of our 
results at /3 = 2.0,2.25,2.5 and 2.75 on 16 3 x 16 lat- 
tice. Some finite-size consistency checks were done on a 
12 3 x 12 lattice. Even though this lattice is relatively 
coarse, we can safely conclude that there are no large 
discretization errors after the implementation of the im- 
proved gauge and fermion actions. 

The glueball correlation function for the ++ channel 
as a function of t at f3 — 2.5 is shown in Fig. ^ The cor- 
relation function exhibits the expected exponential be- 
haviour with respect to Euclidean time t for small values 
of r. Concerning the determination of the tensor glueball, 
we observe that the signal for the E ++ channel is better 
than that for the T 2 ++ channel, and errors on the effec- 
tive mass are reasonable. For this reason, we choose to 
use E of the tensor at finite lattice spacing. At f3 = 2.0, 
the signal does not persist long enough to demonstrate 
convergence to the asymptotic value; we present results 
for time separations t = 1 and 2. We found that Eqs. 
© and (jTU|l yield consistent results in all cases analyzed 
here. The best estimates for the masses from our long 
runs are collected in Table |U 

We extracted the wave functions at time-separations 
t = 1. Little difference between the eigenvectors of C(t) 
at t = 1 and 2 was found. This suggests that there is no 
mixing between states of distinct masses. Typical plots of 
the wave functions, normalized to be unity at the origin, 
for the scalar and tensor glueballs at j3 = 2.25, 2.5 and 
2.75 are shown in Figs. |2Jand|31 respectively. As a guide 
to the eye, Monte Carlo points corresponding to of the 



FIG. 1: Correlation function for the scalar glueball at /3 = 2.5 
on a 16 4 lattice 



same value of are connected with straight lines. The 
scalar wave function exhibits the expected behaviour and 
remains positive for all the values of (3 analyzed here. 
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FIG. 2: SU(3) scalar glueball wave functions computed on a 
16 4 lattice for various values of /3. 

We find that the tensor wave function shows the ex- 
pected flatness and is much more extended than the 
scalar one as (3 increases. This would imply that ten- 
sor is more sensitive to finite-size effects, which is very 
visible in the distortion of the wave function for large r 
at j3 — 2.75. Naively, we expect that the spatial size 
at which we begin to encounter large finite-size effects 
is related to the size of the glueball. From optimization 
analysis at t = 1 /0, we found that the 4x4 and 5x5 loops 
have better overlap with the glueball. This suggests that 
the glueballs have a size of ~ 4a — 5a. The smaller loops 
exhibit a weak signal and a slow convergence with -/V smr . 
This is consistent with the findings of APE in the case 
that a 1 x 1 loop is used as a template . 

The size of the glueball is actually a nontrivial quan- 
tity. Although the charge radius of the glueball can be 
formally defined, its electric charges are carried by the 




FIG. 3: Wave functions of the tensor glueball computed on 
a 16 4 lattice with j3 = 2.25, 2.5 and 2.75. 

quarks and antiquarks, which play only secondary roles 
in describing the glueball state in the idealized limit, 
since the glueball does not contain any valence content 
of quarks (antiquarks). In principle, one could use Eq. 
© to extract the effective radius, but because of the dis- 
tortion of the wave function at large r the results would 
depend strongly on the fit analysis of the above form. 
Even after the ground-state enhancement, the complete 
elimination of all the contributions of the excited states is 
impossible, especially near r ~ 0. It follows that Eq. © 
holds only on a limited interval, which does not include 
the vicinity of r ~ 0. Hence, for the accurate measure- 
ment of the effective radius, we need to find an appro- 
priate fit range. To this end, we examine a plot of the 
effective-radius and the ratio 



for a given ip(r)/i/)(r + 1) at each fixed r. In Fig. 01 we 
plot the effective radius ro as a function of r associated 
with Figs. and El at j3 = 2.75. Figure H shows a non- 
trivial r dependence in the neighborhood of r ~ in the 
effective radius of the tensor state. However, at large r, 
we see that there appears a region where ro takes almost 
a constant value. Thus we conjecture that ro consists of 
nearly a single spectral component, and hence can prop- 
erly represent the effective radius. Estimates of the sizes, 
in lattice units, at various values of /3 values are listed in 
Table HU 

In order to ascertain the finite-size effect on our mea- 
surements, we performed extra simulations for = 2.5 
and 2.75, using a lattice of spatial extent L = 12. The 
results from this spatial extension for glueball masses in 
terms of lattice units are given in Tables HI and ITT1 respec- 
tively. Note that the results for the smaller volume differ 
very little from those for 16 3 lattice, indicating that sys- 
tematic errors in these results attributable to the finite 
volume are negligible. We also find that our estimates 
for the size of glueballs are consistent with the assump- 



FIG. 4: The effective radius for the scalar and tensor glueball 
states at P = 2.75. The dashed horizontal lines indicate the 
plateau values. 

tion that there is no large finite volume dependence at 
f3 = 2.75. Given that we use the larger lattice size to ex- 
tract our results we can disregard any contamination of 
the glueball signal by a torelon-antitorelon pair, since no 
mass reduction of sufficient magnitude was found as the 
lattice volume was reduced (see Table QJ. This suggests 
that none of our states could be interpreted as a torelon 
pair. 

TABLE I: Scalar and tensor glueball energy estimates in 
lattice units for the spatial extensions L = 12 and L — 16. 
Also shown are our estimates of the rho mass at a physical 
point on the L = 16 lattice. 







Mass 






am ++ 


am 2 ++ 


am p (K c ) 


p/L 


12 16 


12 16 


16 


2.0 


0.801(5) 


1.42(2) 


0.563(3) 


2.25 


0.58(1) 


1.07(2) 


0.391(4) 


2.5 


0.52(2) 0.527(3) 


0.82(2) 0.824(16) 


0.288(4) 


2.75 


0.43(2) 0.443(4) 


0.64(2) 0.659(14) 


0.227(4) 


TABLE II: Effective radii 


of scalar and tensor glueballs in 


lattice units for the two spatial extensions L — 12 and L — 16. 






Size 






0++ / 

T-o /a 


2++ / 

r Q /a 


P/L 


12 16 


12 


16 


2.0 


1.37(7) 


3.08(14) 


2.25 


2.9(1) 




5.5(5) 


2.5 


4.1(1) 4.23(4) 


7.6(8) 


7.5(6) 


2.75 


4.8(2) 4.77(6) 


10.0(1.3) 9.95(1.24) 



To convert scalar and tensor glueball sizes and masses 
to physical units and extrapolate to the continuum limit, 
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we first need an estimate of the lattice spacing. One natu- 
ral choice for this conversion factor is the rho mass, m p a. 
Estimating the lattice spacing can be done by extrapo- 
lating the p mass to the physical quark mass 2 . Using a 
tadpole-improved clover fermion action we obtain the es- 
timate 1/k c = 6.362(3) by extrapolating the square pion 
masses from the largest five k values at (3 = 2.5. We find 
that a linear fit of 1/k to m 2 in 1/k works well for the rest 
of the coupling values analyzed here. Linearly extrapo- 
lating the p mass to n c , we find am p (K c ) — 0.288(4) or 
a -1 = 2.66(5) GeV, where the error is a jackknife esti- 
mate. The lattice spacings at other values of (3 are listed 
in Table ITTT1 

Another type of error that might effect the simulation 
results comes from the scaling violation for our action. 
Expecting that the dominant part of the scaling viola- 
tion errors from the gauge and light quark sectors are 
largely eliminated by the tadpole improvement, we ex- 
trapolated the results obtained with finite a to the con- 
tinuum limit, a — > 0. In practice, it is often difficult to 
quantify the magnitude of the systematic errors arising 
from this origin. Here, we adopt an a 2 - linear extrapo- 
lation for the continuum limit, because the leading-order 
scaling violation is always 0(a 2 AQ(j-£ ) m q ). We also per- 
formed an a-linear extrapolation to estimate systematic 
errors. We use the results obtained with the three finest 
lattice spacings for the continuum extrapolation, exclud- 
ing the results obtained at (3 = 2.0, which appear to have 
large discretization errors as expected from the naive or- 
der estimate. 



TABLE III: The lightest SU(3) glueball energy and size es- 
timates in terms of the p mass. 



p 


a(fm) 


m ++/m p 


m 2++ /rn p 


0+ + 

*o rn p 


2 + + 

r m p 


2.0 


0.1444(6) 


1.422(8) 


2.53(4) 


0.77(5) 


1.73(4) 


2.25 


0.100(2) 


1.50(2) 


2.76(6) 


1.15(4) 


2.17(9) 


2.5 


0.073(3) 


1.83(3) 


2.86(4) 


1.18(3) 


2.2(1) 


2.75 


0.058(6) 


1.95(4) 


2.87(5) 


1.12(5) 


2.26(14) 



Performing such extrapolations for mass and size we 
adopt the choice that yields the smoothest scaling be- 
haviour for the final value, and use others to estimate 
the systematic errors. 

Figure displays our results for the scalar and tensor 
glueball radii in terms of m p . It is seen that linear ex- 
trapolations in a 2 yield better fits and the value in the 
continuum limit of 1.14 ± 0.08 and 2.27 ± 0.03 for the 
scalar and tensor states, respectively. 

Note that the product r^ 2 m p varies only slightly over 
the fitting range. The three non-zero lattice spacing 




FIG. 5: Continuum limit extrapolation of glue radii in terms 
of p mass. The dashed curves are the best fits to the simula- 
tion results over the range 0.0033 < a 2 < 0.01. 



values of the product are within 0.02 - 0.04 and 0.01 
- 0.09 standard deviations of the extrapolated zero lat- 
tice spacing result for the scalar and tensor state, re- 
spectively. This allows unambiguous and accurate con- 
tinuum extrapolations. Figure plots the dimensionless 
ratio rao^/m/) as a function of the square of the lattice 
spacing. Again, we expect a linear fit in a 2 to provide 
the most reliable extrapolation to the a — > limit. Lin- 
ear extrapolations to the continuum limit yield a scalar 
estimate of 2.19 ± 0.06 and a tensor result of 2.95 ± 0.12. 
In contrast to the tensor, the scalar glueball mass shows 
significant finite-spacing errors. The continuum limit re- 
sults obtained, in terms of the rho mass, are summarized 
in the Table E| 




FIG. 6: Continuum limit extrapolation of the glue energy 
estimates in terms of the p mass. The curves represent linear 
fits to the data over the range 0.0033 < a 2 < 0.01. 



2 Vaccarino and Weingarten noted that extrapolating to zero To obtain the masses and the radii in units of MeV 
lattice spacing using [A (i ^/g a 0) ] should give results nearly equiv- and fm, we used the experimental value 768 MeV for 
alent to those found using m p (a)a. m p . This yields radii of 0.29 ±0.02 and 0.58 ±0.01 fm for 
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TABLE IV: Continuum limit predictions for the scalar and 
tensor glueball masses and sizes and their conversion to MeV 
and fm, respectively. 



m ++/rn p 

m 2++ /m p 

o++ 
r m p 

2+ + 

r m p 

m 0++ 
m 2 ++ 
.o° ++ 



2.19 ±0.06 
2.95 ±0.12 
1.14 ±0.08 
2.27 ±0.03 
1680 ± 46 MeV 
2265 ± 92 MeV 
0.29 ± 0.02 fm 
0.58 ±0.007 fm 



actions. Instead of fixing a gauge or a path for gluons, the 
correlation functions were computed from gauge invari- 
ant closed-loop paths. This approach has an advantage 
over the method used in Ref which can involve a 
mixing of the glueball states with torelons, and detan- 
gling of the glueball and torelon was done by hand. The 
iterative smearing procedure discussed here can be used 
to obtain a rough estimate of the physical glueball size 
in terms of a Gaussian extension. Although easy to im- 
plement, it does not give detailed information about the 
glueball wave function. A better approach may be to con- 
struct an elaborate smearing procedure that maps onto 
the wave function in one step (i.e., a variational method 
combined with fuzzing techniques). 



the scalar and tensor glueballs, respectively. Our results 
show that the scalar glueball has a radius roughly equal 
to that of a pion [T^j and the tensor glueball is about 
two times as large as the scalar glueball. We thus find 
that glueball size comparable in those of other hadrons. 
The continuum limit glueball mass results in MeV are 
summarized in Table IIVI Our glueball mass results are 
in good agreement with those of previous calculations us- 
ing the improved gauge action [3j as well as with those 
obtained with the Wilson action Q. It is encouraging 
that the results are, within 10 — 15% errors, independent 
of the lattice action. We find that our glueball size esti- 
mates are consistent with the results reported in Ref. 0] 
obtained by using a Gaussian extension as the character- 
istic size of the glueball. We believe that our physical 
glueball sizes arc more accurate than those obtained in 
Ref. where it was found that the radius of the tensor 
glueball is three times that of a pion and four times that 
of a scalar glueball. The predicted zero lattice spacing 
results are in fact not found by extrapolation to zero lat- 
tice spacing, but are instead obtained from calculations 
with (3 = 2.2 for the glueball effective radius. Assigning 
a physical value to the lattice spacing a by setting the 
string tension to 420 MeV, the physical radius is calcu- 
lated with zero uncertainty, with no accurate representa- 
tion of the effect of the absence of extrapolation. 



IV. CONCLUSION 

To conclude, we have calculated the SU(3) glueball 
wave functions using improved gauge and fermion lattice 



Finally, we note that our results at the smallest lat- 
tice spacings seem to scale. This allows for very accurate 
extrapolations to the continuum limit. It is found that 
the physical size of the tensor glueball is about two times 
as large as that of the scalar glueball. This is consis- 
tent with the large finite-size effect for the tensor (E ++ ) 
glueball mass found in previous lattice calculations. A 
great deal of care should be taken in making direct com- 
parisons with experiment, since these values ignore the 
effects of light quarks and mixings with nearby conven- 
tional mesons. We intend to include the effect of dynam- 
ical quarks in a future study. 
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